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ABSTRACT 

Gluon bremsstrahlung in scattering events with high transverse momentum jets is expected 
to increase markedly with the hardness (J2Et) of the primary event. Within perturbative 
QCD we estimate a probability of order unity to see additional minijets with Et ~ 15 GeV 
in "dijet" events with Et > 400 GeV. The veto of such minijets is a promising background 
rejection tool for the Higgs search at the LHC. 
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The study of detailed properties of high transverse energy (Et) jet events in pp collisions 
at the Tevatron has become an important testing ground for perturbative QCD (pQCD)[|TJ. 
At the same time the investigation of hadronization patterns in dijet events with widely 
separated jets has lead to the discovery of rapidity gaps at the Tevatron |2| and has demon- 
strated that the color structure of the hard scattering event has dramatic consequences for 
the angular distribution of produced hadrons at hadron colliders |3], f|, || . More precisely it 
is the t-channel exchange of color singlet quanta which can lead to rapidity gaps. One of the 
most important members of this class of processes will be weak boson scattering at hadron 
supercolliders, i.e. the electroweak process qq — > qqVV, and, indeed, a rapidity gap trigger 
had been suggested as a promising Higgs search tool at the SSC 

A modified rapidity gap trigger, in terms of minijets of Et > 15 — 20 GeV instead of 
soft hadrons, was recently suggested for the study of weak boson scattering at the LHC ||. 
Due to the high luminosity of the LHC, which leads to multiple pp interactions in a single 
bunch crossing, the search for rapidity regions without soft hadrons is rendered impractical. 
Minijets of sufficiently high Et, however, are unlikely to be produced by multiple interactions. 
On the theoretical side rapidity gaps are not well described by pQCD which makes predictions 
about jets rather than individual hadrons. Minijets, on the other hand, can be described by 
pQCD if the perturbative expansion is taken far enough. 

A key component of the minijet veto is the observation that the emission of semihard 
partons, soft at the level of the hard scattering process but still leading to distinct minijets, 
becomes quite probable in the very hard QCD background processes to weak boson scatter- 
ing. A tree level calculation for the QCD W pair background reveals that the production 
cross sections for WWj and WWjj events become equal when the minimal Et for which 
partons are identified as jets is taken to be of order 30 — 40 GeV ||. This leads us to 
expect the emission of multiple minijets with Et > 20 GeV in these background processes. 
The -Ey-scale can be understood qualitatively by noting that the invariant mass of the two 
W's must be in the 600-1000 GeV range to become a background for a heavy Higgs boson 
and the presence of an additional forward tagging jet increases the required cm. energy 
y/I into the TeV range. The emission of additional gluons above a transverse momentum 

■ 2 

i?™ m is suppressed by a factor f s « a s In (s/E T ^ in ) and this factor approaches unity for 
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E™ m = 0(30 GeV). Multiple minijet emission at such large E^s should be observable even 
in a high luminosity environment and it is the veto of these minijets which should substan- 
tially reduce the backgrounds with little effect on the signal rates (due the different color 
structure and QCD scales of the latter). 

The approximate equality of n- and n + 1-parton cross sections clearly indicates that 
fixed order pQCD is no longer reliable in this range. One would like to have experimental 
data on minijet multiplicities, veto probabilities, angular distributions and the like and to 
then gauge the pQCD calculations for new physics backgrounds in the light of these direct 
observations. In this letter we suggest that such information can already be obtained in the 
hardest QCD events available now, namely in hard "dijet" events at the Tevatron. 

In the following we consider pQCD predictions for dijet production at the Tevatron 
in next-to- leading order (NLO), using the NLO Monte Carlo program JETRAD of Giele, 
Glover and Kosower |7]]. MRSD'_ distribution functions || are used and the strong coupling 
constant a s (fi) is evaluated at 2-loop order with ot s (mz) = 0.111. Both the renormalization 
and the factorization scales are chosen as n — J2Et/A, where J2 Et is the scalar sum of 
the parton transverse momenta. The program allows us to make accurate estimates for 
the two-jet inclusive cross section (o^mci) at the 1-loop level. In addition, three parton 
final states are simulated at tree level which allows the study of soft minijet activity at 
leading order. For the three parton cross section (03) the usual scale ambiguity of tree level 
calculations arises, which in the present case is exacerbated by the very different scales of 
the hard event (J2Et = O (500GeV)) as compared to the transverse momentum of the 
softest parton (Et 3 = O(20 GeV)). Truly we are dealing with a 2-scale problem and a 
rough estimate of the ambiguities is obtained by comparing <7 3 (/i = Y^E^/A) with a 3 (/i = 
J2 Et /4) ■ a s (Et 3 )/a s (J2 Et /4) i.e. by effectively using the softest Et as the scale for emission 
of the additional parton. We will refer to these two choices as ^ = J^Et/A and ^3 = Et 3 , 
respectively. 

We cluster partons into jets using a DO style cone algorithm ||. If two partons are to be 
clustered into a jet we form the jet momentum as 

Pjct = Pi + P2 , (1) 
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E ^et = Et 1 + Et 2 ■ ( 2 ) 

Each pair of partons is clustered into a temporary jet. Whenever both partons are within a 
distance 

Ai?(jet, parton) = yj (r) ict - r] par ton) 2 + (^jet ~ ^parton) 2 < 0.7 (3) 

of the temporary jet, then these partons are clustered into an actual jet. Since we are 
dealing with a maximum of three partons in the final state, at most one such recombination 
of massless partons can occur in our program. 

Final (clusters of) partons are then counted as jets if their transverse energy Et and 
pseudorapidity 77 satisfy the conditions 

E T >E T lin , M<3.5, (4) 

with, typically, Ej? n = 15 GeV. We define the J2Et of the event as the sum of the Et's of 
all observed jets. 

In Fig. [l| we show the fraction of 2 jet inclusive events that contain a third jet, for values 
E™ m = 15, 20, and 30 GeV. The 3 jet fraction / 3 = o^/o^inei is plotted as a function of 
the Et of the event and for the two choices of the effective scale governing soft minijet 
emission. Below ~ 200 GeV one observes a marked rise of the fraction of three jet events 
with Et which then flattens out around Et = 300 GeV and slowly decreases for larger 
values. This pattern can be understood by considering the phase space, both in Et and 
pseudorapidity, which is available to the third parton. By definition, this is the parton with 
the smallest Et and hence its allowed range is E™ m < Et < J2Et/3. Integration of the 
partonic cross section da ~ cIEt/Et over this range leads to a logarithmic increase of the 3 
jet fraction and explains the rise of u^/ <72,\nc\ for small J2Et- 

For 2-jet events and, approximately, for 3-jet events with a soft central jet the rapidity 
difference of the two highest Et jets and the cm. energy are related by 

v / I = ^E T -cosh^^ . (5) 

Due to the rapid decrease of parton luminosity with rising s only small rapidity differences 
rji — 772 can be reached at high J2 Et- This effect is demonstrated in Fig. ^: the rapidity 
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Figure 1: Fraction of 2-jet inclusive events that contain a third jet within the acceptance 
requirements of Eqs. The 3-jet cross section is integrated above a minimal Et and 
above minimal transverse energies Et 3 of the softest jet of 15 GeV (solid line), 20 GeV 
(dashed line) and 30 GeV (dotted line). The renormalization scale is set to fi = J2Et/4= in 
part a). Changing the effective coupling constant for soft parton emission in the 3-jet cross 
section to q: s (Et 3 ) yields the results shown in part b). 

distribution in the cm. frame, da/d\r]i — 772 1 becomes narrower with increasing J^Et- The 
emission of a third jet well outside the rapidity interval set by the two hardest jets will 
markedly increase the s of the event and hence will suffer from a lower parton luminosity. 
Thus additional soft jets are predominantly produced in the rapidity range set by the two 
hard jets and therefore the decrease of the average \r]i — r/ 2 | with increasing reduces 
the phase space for emission of additional jets. This explains the decreasing fraction of three 
jet events at large Et which was observed in Fig. [I]. 

These phase space effects are very process specific. In weak boson scattering events at the 
LHC, for example, (and backgrounds to these events) the available s is largely determined by 
the forward and backward tagging jets and hence the kinematically favored rapidity range 
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Figure 2: Normalized rapidity distribution I/02 • da2/d\rji — 772 1 for 2-jet inclusive events at 
the Tevatron. Shown is the narrowing of the rapidity range \rji — r] 2 \ with increasing hardness 
E T ) of the event. 

for minijet emission is much larger than in dijet events at the Tevatron. For a comparison 
of different processes one would like to eliminate these purely kinematical effects. For jet 
production at the Tevatron our Monte Carlo studies show that, to good precision, the prob- 
ability for a 3rd jet to be produced in the 77 range between the two hardest jets rises linearly 
with the r]-cf) phase space which is available for that third jet, that is 

^(V3 e [771,772]) 



A 



CYincl 



lego 



27r|^i - 772 1 - A c 



(6) 



where A cone is the part of the lego-plot region between 7/1 and 7/2 which is covered by the two 
cones of radius 0.7 about the two hardest jets. We therefore divide out the available phase 
space A\ ego and consider the probability density for emission of a third jet in the lego-plot. 
In Fig. ^ we show the average density 

1 d 2 a 3 
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for emission of the third, lowest Ej<, jet in the interval [771, 772] set by the two hardest jets. The 
probability density as a function of the minimal Ex of the event is shown in Fig. Rk, while 
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Figure 3: Lego-plot density pi ego of soft jets between the two hardest jets. Results are shown 
a) as a function of the minimal Et and for three values of E™ m and b) as a function of 
E™ m with Et as a parameter. The effective coupling constant for soft parton emission in 
the 3-jet cross section has been set to a s (-&r 3 )- On the r.h.s. the scale is multiplied by a 
factor 2-7T to show the approximate density per unit of pseudorapidity. 

Fig. |3|b depicts the same probability density as a function of E^ m , the minimal transverse 
momentum of the minijet. One now observes a monotonic increase with the hardness of the 
overall event, Ex- It is this increase which is characteristic of the effect of enhanced minijet 
emission and which allows to distinguish it from backgrounds due to the underlying event. 
Factorization of the hard event implies that the background from the underlying event, that 
is minijet emission caused by the interactions between the colored remnants of the original 
pp pair, does not grow as the X) Ej- of the hard event increases. The minijet background 
from the underlying event may even decrease at very high Et because less longitudinal 
momentum is available for the interactions between the pp remnants. 

For very hard events, Et > 400 GeV say, the minijet emission probability reaches 
remarkably large values, of order 0.3 per unit of rapidity for minijets of Et > 15 GeV. 
These large probability densities reflect the approximate saturation of the 2-jet inclusive 
cross section with 3-jet events which is evident in Fig. |I[ Clearly, for hard jet events at 
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the Tevatron with J2 Et ~ 400 GeV the calculation of jet multiplicities within fixed order 
pQCD is breaking down for jet thresholds as large as i?™ m = 30 GeV. One should expect large 
probabilities for 4 or 5 jet events as well, but a reliable prediction of these multijet fractions 
goes beyond the perturbative tools at hand. A leading logarithm parton shower calculation 
does allow to simulate such events. However, the scale of the logarithms remains uncertain 
and the multijet rates are again subject to large normalization uncertainties. Experimental 
input is needed in addition to theoretical considerations. It should be noted that these events 
are not rare, indeed ct(J2Et > 400 GeV) ~ 200 pb. Thus we expect a large sample of high 
J2 Et events at the Tevatron which exhibit multiple emission of additional minijets in the 
10 — 30 GeV Et range and which can be used to study the characteristics of minijet emission. 
Specific questions which should be addressed experimentally are 

1. What is the average multiplicity of minijets and how does it change with the YE? of 
the event? 

2. How are minijet multiplicities distributed? The NLO Monte Carlo only allows to 
generate events with and 1 minijet. 

3. Down to which E™ m can one distinguish minijet activity arising in the hard scattering 
event from minijets in the underlying event? 

4. What is the probability for zero minijet emission in a hard scattering event? How does 
this probability change with Y E^l How does it depend on the minimal transverse 
energy E™ m which is required for a cluster of hadrons to be identified as a jet? 

5. Does the effective coupling, which governs the emission of minijets, depend on the hard 
scattering scale, the soft scale, or some combination of the two? 

The results shown in Fig. [I] can only be considered crude estimates, as can be seen 
from the differing predictions when changing the scale of soft emission from Et/4 to Et 3 - 
Theoretically this question of scale can be resolved by going to the next order in pQCD, 
in this case calculating 3 jet production at NLO. Experimental information on this scale 
is obtained by comparing the rise of the lego-plot density as a function of J2 Et with the 
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/z 3 = E T3 prediction in Fig. |a. Considering the large uncertainties of our pQCD calculation 
it is well possible, however, that a relatively large correction factor F is needed to relate the 
experimental minijet density with the pQCD expectation, 

F r 2n 1 d 2 a 3 
Minijet density = Fp lego = / dr) 3 dfo - — — . (8) 

o"2,inci Jm Jo A lcgo di] 3 d(p 3 

Measuring the dependence of this correction factor on the kinematics of the hard event (J2 Ft 
and 1 771 — 772!) would provide a useful guide for theoretical improvements, which would then 
help to more reliably predict minijet emission probabilities in hard scattering events at the 
LHC. 
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